Dark entangled steady states of Rydberg interacting atoms 
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We propose a scheme for rapid generation of high fidelity steady state entanglement between a pair 
of atoms. A two-photon excitation process towards long-lived Rydberg states with finite pairwise 
interaction, a dark state interference effect in the individual atoms, and spontaneous emission from 
their short-lived excited states lead to rapid, dissipative formation of an entangled steady state. We 
show that for a wide range of physical parameters, this entangled state is formed on a time scale 
given by the strengths of coherent Raman and Rabi fields applied to the atoms, while it is only 
weakly dependent on the Rydberg interaction strength. 

PACS numbers: 42.50.Gy, 42.50.Dv, 03.67.-a, 42.50.-p 



The strong blockade interaction between Rydberg ex- 
cited atoms open many possibilities to explore neutral 
atoms for quantum computing and for the study of a va- 
riety of complex many-body and light-matter problems 
[l|. The first proposal by Jaksch et al.^ to use Ryd- 
berg blockade to implement a fast two-qubit controUed- 
NOT (CNOT) gate has been followed by a variety of 
schemes for fast quantum gates with atomic ensembles 
[S-Q, entangled state preparation [7|, quantum algo- 
rithms [8|, |9| , quantum simulators |lO| , and efficient quan- 
tum repeaters |ll|. Fidelities of around 0.9 for perform- 
ing a CNOT gate (probability truth table) and 0.7 for 
generation of entanglement between two atoms using the 
Rydberg blockade interaction have been reported in jl2| - 
Il4l | . These fidelities are mainly limited by the finite mag- 
nitude of the blockade interaction with respect to the ex- 
citing lasers and errors due to spontaneous emission from 
the intermediate state used in the two-photon coupling to 
the Rydberg state [l5|. To reduce these errors one must 
excite very high lying Rydberg states with large blockade 
interactions and one must apply excitation fields with a 
large intermediate state detuning. Dissipation provides 
an attractive supplement to unitary interactions for the 
preparation of few-atom entangled states and can even 
be tailored to implement universal quantum computa- 
id llg . 11?! . This has drawn further attention to the 
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use of dissipation as an active ingredient in quantum in- 
formation processing. In this Letter, we show that decay 
by spontaneous emission of light can be used as a key 
resource to generate pairwise entangled steady states of 
atoms with interacting Rydberg states. We shall show 
that our scheme works for even moderate interactions 
between the excited atoms and that the convergence to 
the entangled steady state is rapid enough to yield ro- 
bustness against realistic noise and loss mechanisms. 

The physical setup for our dissipative generation of 
entangled states involves two atoms with two ground hy- 
perfine states, labeled |0) and |1), and a Rydberg state 
|r), which can be excited via the intermediate state \p), 



FIG. 1: Schematic representation of two atoms driven by laser 
and Raman fields {Q.i,Q.2,u)) and coupled via Rydberg state 
interaction Vrv- The atoms are described as four- level systems 
comprised of two ground levels |1) and |0), coupled weakly to 
each other by a Raman process, while the ground level |i) is 
strongly coupled to the Rydberg level \r) by a resonant two- 
photon process via the intermediate short-lived excited level 



see Fig. 1. The atoms are trapped within a distance 
of few tens of /im, such that they experience a non-zero 
energy shift 14-r, when both atoms occupy the Rydberg 
state \r). We apply resonant excitation from the state 
|1) to the optically excited state \p) and from \p) to the 
Rydberg level \r) with Rabi frequencies fii and il2, re- 
spectively, and we drive the transition between states |0) 
and |1), by a resonant Raman process with strength w. 
A detuning A may be applied with respect to the inter- 
mediate level |p), but in the following we will set A = 0. 
The Rydberg and optically excited states decay by spon- 
taneous emission of radiation. We assume that the Ryd- 
berg state lifetime is much longer than the optical state 
lifetime 7^ << 7p. 

To briefly describe the main idea behind our entangle- 



ment mechanism, it is useful to consider first the states 
of a single atom in the absence of the Raman coupling 
field. The laser fields give rise to the electromagneti- 
cally induced transparency (EIT) phenomenon [18| asso- 
ciated with the coherent trapping of atomic population 
in the dark eigenstate \D) = ^\p,2\l) — ^i\r)], where 
Q, = -y/fi^+Tif . The atom may also reside in the uncou- 
pled ground state |0)), and if we apply a weak resonant 
Raman coupling between states |0) and |1) with coupling 
strength w [lo < r2i,J72), we effectively couple the pair 
of states |0) and \D) with the coupling strength uji^2/^ 
given by the state |1) component of \D). 
The dark singlet combination, 



\DS)^^i\DO)-\OD)), 



(1) 



of two atoms is invariant under the application of the 
strong Rabi and the weak Raman coupling. The triplet 
spaceof states (|00), {\D0) + \0D))/V2, |D£))) are also 
dark states, i.e., they do not couple to the short lived 
excited atomic state, but the Raman coupling causes ro- 
tations among them. In the absence of interatomic inter- 
action, any initial product state of the atoms, will on the 
time scale of a few 7"^ evolve into an incoherent mixture 
of the dark singlet and triplet states with no entangle- 
ment. 

Since the singlet state contains only one atom in the 
state \D) and hence only single Rydberg state compo- 
nents, it is immune to the interaction term Vrr, while 
the triplet component \DD), due to the interaction, is 
energetically shifted away from the ground state man- 
ifold and seizes to be an eigenstate of the interaction. 
This destroys the EIT mechanism [19I-I2H and the |11) 
component of \DD) now couples to the rapidly decaying 
states. As a consequence of this coupling, the state |11), 
and the other two triplet states, which are coupled by the 
Raman field, acquire finite lifetime and undergo excita- 
tion and decay by spontaneous emission until the atomic 
population accumulates in the dark singlet state ([T]). Ide- 
ally, when the blockade strength becomes infinite, all the 
population will rapidly return into states with maximally 
one atom in the state |D). The presence of finite block- 
ade has the same effect, but may require a longer time 
(^ 'iil.'^ /nIVrr) to remove also the population from the 
doubly excited long lived Rydberg states. When the asso- 
ciated rate is larger than the Raman field strength w, the 
system does not distinguish between finite and infinite 
blockade strength, and contradictory to the condition for 
unitarily generated entanglement, which strictly requires 
Vrr- 3> rj, the much weaker condition 
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will ensure the formation of a high fidelity entangled 
state. This entangled steady state is obtained from any 
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FIG. 2: (Color online) The fidelity F{t) = {DS\p{t)\DS) (red 
dashed line), the purity of the two-atom state P{t) (green dot- 
dashed line) and the population of the effective four-level sub- 
space Peff{t) (blue solid line), are plotted as functions time. 
Shown in the inset is the early time variation of the same 
quantities. An analytical estimate for the fidelity Feffit) 
(black dotted) obtained from the effective four-level subspace 
is also shown in the figure. The parameters chosen for the 
calculation are ni/27r = 20 MHz,r22 = 2r2i, tj/27r = 250kHz, 
7i/27r = 72/27r = 3.03 MHz, Kr = ^i and 7r/27r = IkHz. 



initial state of the atoms, assuming neither a very large 
nor a particularly precisely defined value of the interac- 
tion. 

For a quantitative analysis of the dissipative prepara- 
tion scheme, we shall solve the two-atom master equation 
numerically and we shall investigate and interpret the de- 
tailed dependence of our results on the physical param- 
eters of the problem. The total Hamiltonian describing 
the interacting Rydberg atoms is 

H = Hi®I + I®H2 + Vrr\rr){rr\, (3) 

with {j = 1,2) single atom Hamiltonian operators, 

-^^\\)^^{p\+^2\p),,{r\+h.c, (4) 



i^=-|l>,,< 



where we assume a resonant coupling to take full advan- 
tage of the rapidly decaying intermediate level \p) . 

The evolution of the system due to spontaneous emis- 
sion is described by the master equation 



9,p = z[He//,p]+E4'^ PC^^^ Heff=H-^Y.^i^^ 4^ 
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where C^^' are Lindblad operators, which describe the 

(2) decay processes of the j — th atom by spontaneous emis- 
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sion of light, Cq 

ilp ~ 70 + 7i). Rydberg state decay is described by 

similar terms, e.g., Cr — \/lr\p) jj{A^ but since the Ry- 
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FIG. 3: (Color online) The dependence of tlie fidelity F(t) 
on time is shown for four different values of Yrrj^ keeping 7p 
fixed. In the inset is shown the variation of F{t) for different 
values of the decay rate 7p for perfect blockade condition. For 
the above simulation the other fixed parameters ni/27r = 20 
MHz,!:!2 = 2^1, a;/27r = 125kHz, 7i/27r = 72/27r = 3.03 MHz 
and 7r/27r = IkHz. 



dberg state is long lived, we defer our discussion of its 
consequences for the scheme to the end of the Letter. 

In Fig. 2 we have plotted the fidelity F[t) — 
{DS\p{t)\DS), and the purity P{t) =Tr(p(i)2) of the 
state determined by solution of the two-atom master 
equation, starting from the product state |11) (the results 
are similar for any other initial state). The results con- 
firm the evolution described above, with an initial rapid 
evolution towards an entangled state fidelity of approx- 
imately 0.25 and a comparable low value of the purity, 
compatible with the atoms incoherently populating prod- 
uct states of |0) and |D) (see the inset for a magnified 
view of the initial dynamics) . 

To explain the results in Fig. 2 we shall describe the 
long time behavior within the effective four-level sub- 
space spanned by the basis vectors ITq) = |11), l^^i) = 
^[\D0) + \0D)], \T,) = 100), IT3) ^ \DS) = j^[\DO) 

\QD)]. The effective unitary evolution can be found by 
projecting the total Hamiltonian (|3]) on to this subspace, 
which gives 

n^Cj {\To){T,\ + |ri)(r2| + h.c.) , u;^V2u: (^^ .(6) 

In addition to the pure unitary evolution there is a strong 
coupling of the level |To) to the high energy spectrum. 
The rapid decay from these high energy states effectively 
appear as if the population in |To) is decaying back into 
all four states of the subspace, and the effective decay 
rate F will be determined by the Rabi fields, 7^, and ¥„■ 
The new Lindblad operators (assuming perfect blockade) 



C, = ]Jj\Tj}{To\ 



'^ {-^) 



(7) 



where j = 0, 1, 2, 3. As jTs) is a zero eigenvalue 
state of the Hamiltonian and a dark state with respect 
to the dissipative dynamics, it accumulates the popula- 
tion on a time scale determined by w and F. In Fig. 2 
one can see an excellent agreement between the exact 
numerical simulations and the dynamics predicted by 
the effective four-level dynamics described here. For 
the parameters used in the plot, 7p/27r = 6.07MIIz, 
ri2/27r = 2ni/27r = 40MHz and uj/2tt = 0.25MHz, the 
effective decay rate and frequency are F/27r — 1.89MIIz, 
and a)/27r = 320kIIz. To further confirm the effective four 
level dynamics we have plotted the population of the ef- 
fective four level subspace Pe// (i) obtained from the full 
master equation analysis. One can see that all the pop- 
ulation resides in the four level subspace for t > 1/F. 
Hereafter, on a longer time scale the system reaches a 
steady state with a high overlap with the dark singlet 
state \DS), ^. 

As long as the coupling and interaction parameters ful- 
fil fli, ^2, 7, Vrr ^ uj and $^2 > ^1, Fig. 3 reveals a very 
weak dependence of the fidelity F{t) on the dissipative 
atomic decay rate 7^ and the interaction strength Vrr (for 
more information see |22| ) . From the same figure one can 
see that for Kr/f^ > 1/4, and jp/ui > 20 {jp/il. > 0.1), 
the formation of the entangled state is almost purely de- 
termined by the interplay of dynamics between the Rabi 
and Raman fields. This is surprising since spontaneous 
decay and Rydberg interactions are both essential physi- 
cal components for the formation of the entangled state. 
Only when either of these strengths are quite small in 
comparison to the Rabi fields, the convergence to steady 
state becomes slower. In addition to the deviations from 
these conditions the fidelity also gets reduced by other 
long time effects such as decay of the Rydberg state which 
reduces the Rydberg state decay reduces the steady state 
fidehty -F(oo) sa 1 — 7r/F [2^, and dephasing due to e.g, 
magnetic noise. The bounds on the steady state fidelity 
and the rate of entanglement can be determined from 
the eigenspectrum of He// (6). The maximum fidelity 
achievable is given by the overlap of the state \DS) with 
the eigenstate of H^ff which has an eigenvalue with the 



smallest imaginary part 22|. As the imaginary part of 



the energy determines the decay rate for the correspond- 
ing eigenstate the convergence towards the state \DS) is 
determined by the gap between the two eigenvalues with 
lowest imaginary component [2^ . 

In summary, we have presented a fast and robust 
scheme for dissipative generation of entangled steady 
states of a pair of atoms with interacting Rydberg states. 
The scheme does not demand the interaction to be in the 
blockade regime, and provided an experimentally natural 



hierarchy between the magnitudes of the different cou- 
pling terms, it does not depend on fine tuning of any pa- 
rameters. Note that the role of Vrr is merely to perturb 
the state \DD) so that the triplet states are no longer 
dark states of the system. For realistic parameters used 
in the simulations we see that one could achieve a mini- 
mum of 98% fidelity within few tens of microseconds. 

The present analysis also shows how we regain the EIT 
dark state feature (which is lost in the presence of inter- 
action for the case of three level atoms [l^l) by coupling 
weakly to an additional ground state. The effective two- 
level picture of the four level atoms leading to the iden- 
tification of singlet and triplet-like states of two atoms, 
may find applications when studying similar effects in 
larger number of Rydberg interacting atoms. 

Our scheme works under steady driving conditions, 
and by the incorporation of dissipation, it protects the 
desired entangled state and automatically compensates 
for the harmful effects of decoherence and decay. When 
the applied fields are switched off the atoms can be 
driven adiabatically or by a resonant drive into an en- 
tangled states involving only the ground states |0) and 
|1) ([|10) — |01)]/-\/2), and it may then be used as a tele- 
portation channel to accomplish quantum gates between 
remote qubit atoms in the vicinity of the entangled atoms 
in neutral atom arrays. 

This work was supported by the project MALICIA un- 
der FET-Open grant number 265522, and the lARPA 
MQCO program. 
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SUPPLEMENTARY INFORMATION 

We have shown in the main text that after a nonex- 
ponential transient behavior at short time t < l/7p, the 
fidelity F[t)^ converges towards unity in an exponential 
manner. 



F{t > 7p) = Frna.a - dc-''^^''^^) , 



(8) 



where 1 - d = {DS\p{t)\DS) is the population in \DS) 
after the initial transient time r, je is the rate of for- 
mation of the entangled state, and Fmax is the max- 
imum attainable fidelity for given physical parameters 
(rii,ri2,a;,7p,K-r)- 

Rate of entanglement formation - In Fig. 4 (Left) we 
show the variation of je as a function of w keeping all 
other parameters fixed. The results shown are obtained 
by numerical solution of the master equation (Eq. 5, 
main text), and one clearly sees a predominantly linear 
dependence on small co, while a deviation is seen for larger 
values of oj. The calculations also show that if 7p is in- 
creased proportionally with cj the linear behavior holds 
even for large to. 

In Fig. 4 (Right) we show the variation of •je as a 
function of fii, keeping the other parameters fixed. One 
can see that je depends on the ratio e = 172/^17 while it 
is independent of the absolute magnitudes of the Rabi- 
frequencies fii and Q2- Using these observations and 
assuming values for jp and Vrr in the saturation regime 
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FIG. 4: (Left) The entanglement rate 7b is plotted as a function of the Raman frequency lj, for three different cases (as shown 
in the legend). Other parameters used in the simulation are ili/2-K = lOMHz, il2 = Sfii, Vrr = f^i and a;o/27r = IMHz, 
7o = 6MHz. (Right) The entanglement rate is plotted as a function of the Rabi frequency Sli for two different cases: (i) $12 is 
kept fixed (blue squares) and (ii) $12 is varied keeping r22/r2i fixed (red circles). For this plot, Jp/2tt = IMHz, uj/2ir = O.lMHz. 



(see Fig. 3a of main text), the rate of entanglement for- 
mation can be approximated by 



Bounds on the rate and magnitude of entanglement generated 
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'lh,Vrr), 



(9) 



where the function g{j, Vrr) is a constant of order unity 
for a wide range of parameters of the model. 

Maximum fidelity - We now analyze the asymptotic fi- 
delity achievable with the proposed scheme. In the case 
non-interacting atoms, J7i determines the rate at which 
population accumulates in the effective two-level sub- 
space spanned by the single- atom dark state \D) and 
|0). If $^2 ^ ^ij the ground state component of the 
dark state increases but at the same time the conver- 
gence rate to the entangled steady state decreases. If, 
conversely, fli > 0,2, the dissipative processes suppress 
the sensitivity to the Rydberg interaction Vrr, mediated 
by fl2 on the \p) — ^ \r) transition, and f22 > f^i is a 
quite stringent requirement to achieve the high fidelity 
singlet state. In Fig 5 (Left) we present numerically cal- 
culated values of F{t) and Fmax for different values of 
the Rabi frequencies. As shown in the inset, the singlet 
state fidelity is governed by e, with the optimal choice 
for both the convergence rate and for the steady state 
fidelity given by 1.7 < e < 2.3. 

In Fig. 5 (Right) we show the dependence of Fmax on 
the Raman coupling w and the Rydberg state life time jr- 
Clearly, there exists an optimal regime in the parameter 
space to achieve fidelities close to unity. Though a finite 
Vrr is enough to produce the state \DS), a small value 
of Vrr *^ ^ (see Fig. 3a main text) will slow down the 
entangling process. This would in turn increase the errors 
due the finite lifetime of the Rydberg state |r), leading 
to a decrease in the value of Fmn^. 
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resonant couplings 


the Hamiltonian 


(Eq. 


3 main text) has 


three dark states. 
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fiiwlTpo) — (ri| — aj^)|00)], where Ci are the respective 
normalization constants and \Sij) = ^[|*j) ~ b*)] 
and \Tij) = -75[|jj) + \ji)] are the respective singlet 
and triplet combinations of the basis vectors \ij), and 
hence any linear superpositions of the three states 
also remain dark. In the presence of dissipation 
neither of these remain dark, as each of them has 
a nonzero overlap with the decaying level, |p), and 
hence will gain an imaginary component, allowing 
any population in these states to decay in the long 
time limit. While states \D^) and -75I-D1+-D2) attain 
an imaginary component ^ 7(17i/r2) to their eigen- 
value the other dark state |-D-) ~ -j^\F>i — D2) ~ 

:[\Sro) - \Sw) + 3^/2(|5ip) - |5p,))], where 



C — (ip^); gsts the smallest imaginary component 
^ 2^^. Hence \DJ) ~ \F)S) will remain dark even in 
the presence of dissipation. In fact it is closest to the 
maximally entangled supersinglet 

The bounds on the steady state fidelity and the rate 
of entanglement can be determined from the eigenspec- 
trum of "He// (Eq. 5, main text). The maximum fidelity 
achievable is given by the overlap of the state \DS) with 
the eigenstate of -ffe// which has the eigenvalue with the 
smallest imaginary part. As the imaginary part of the 
energy determines the decay rate for the corresponding 
eigenstate the convergence towards the state \DS) is de- 
termined by the gap between the two eigenvalues with 
lowest imaginary component, assuming perfect blockade 




FIG. 5: (Left) Fidelity is shown as a function of time for three different ratios of the Rabi frequencies. In the inset we have 
shown Fmax as function of Rabi frequencies Sli, ^2. For this plot we have taken uj/2n = 250kHz, 7i/27r — 72/27r = 3.03MHz 
and Vrr/'i'T ~ 20MHz. (Right) The maximum value of fidelity, Fmax obtained in the long time limit is plotted as a function 
of iu and 7^. For this plot we have assumed perfect blockade and ijJo/2tt — IMHz, "fro/2n = lOkHz. The values of other 
parameters used in this plot are similar to those in Figure 4. 



the rate of entanglement 



7-B < ^- 
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and the maximum attainable fidelity is given by 
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